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Abstract 

Using weak coupling methods McLerran and Venugopalan [l| expressed the parton 
distributions in large nuclei as correlation functions of a two dimensional Euclidean 
field theory. The theory has the dimensionful coupling g^^, where /i^ ~ A^/^ is the 
valence quark color charge squared per unit area. We use a lattice regularization 
to investigate these correlation functions both analytically and numerically for the 
simplified case of SU{2) gauge theory. In weak coupling {g^j^L << 5), where L is the 
transverse size of the nucleus, the numerical results agree with the analytic lattice 
weak coupling results. For g^^L » 5, no solutions exist at 0(0"^) where a is the 
lattice spacing. This suggests an ill-defined infrared behavior for the two dimensional 
theory. A recent proposal of McLerran et al. for an analytic solution of the 
classical problem is discussed briefly. 
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1 Introduction 



In Ref. ||^] McLerran and Venugopalan proposed that weak coupling methods can be 
used to compute small x parton distribution functions in large nuclei. They wrote 
down a partition function for wee partons with x « A"^/^ in the presence of external 
sources which are the valence quark charges. The only large component of the valence 
quark current is J"*", which is modelled by 

J!: = 5^'+Paix+,x^)6ix-), (1) 

where pa is the density (per unit area) of valence quark color charges. Their partition 
function is obtained by integrating the QCD partition function, coupled to the above 
static current, over all paS with a Gaussian weight. The variance of this Gaussian 
distribution of valence quark charges, p"^ ~ average valence quark 

color charge squared per unit area, is the only dimensionful parameter in the theory. 
If >> Aq^£,, as(/i^) << 1 and weak coupling methods can be used. This model 
could then be studied as a toy model to understand both the rapid growth of structure 
functions at small x and the eventual saturation of these structure functions as 
dictated by unitarity Note that the model of Ref. is gauge invariant due to 

the Gaussian distribution for the valence quark densities. 

In Ref. |5| the saddle point solution of the partition function in the presence of the 
Gaussian random source was obtained by solving the classical Yang-Mills equations 
Df^F^'^ = . Here is the covariant derivative, F^^'^ the non-Abelian field strength 
tensor and is the current in Eq. ||. It was shown that the classical background 
field that satisfies the Yang-Mills equations has a simple structure. Consequently, 
the classical parton distributions can be expressed as correlation functions of a two 
dimensional Euclidean field theory. This is not too surprising since it is well known 
that at very high energies the longitudinal and transverse coordinates decouple. In- 
deed, it has been proposed recently that the limit of x ^ and color ^ oo \s an 
exactly solvable two dimensional field theory Q. In papers subsequent to Ref. ||5|, 
the problem of quantum fluctuations about the background field f^-Q and that of 
initial conditions in heavy ion collisions were addressed 1^. For a brief review of 
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these results, we refer the reader to Ref. [12]. An excehent introduction to all aspects 



of the low X problem is given in Ref. 

In this paper, we will discuss only the classical solutions of the Yang-Mills equa- 
tions. As we shall see in Section 2, computing the correlation functions requires that 
we solve a stochastic differential equation for each color charge configuration. Since 
the equations are highly non-linear, no analytic solutions were found. However, it 
was claimed in Ref. Q that the parton distributions have the Weizsacker-Williams 
behavior in the weak coupling region asfJ- << h: dN/dxd?kt oc 1/xk^. It was conjec- 
tured that the solution of the stochastic differential equations in the strong coupling 
region of Aqcd « osfJ- « /U would reveal that the classical gluons generate a 
screening mass niscreen ~ osA*- If there is such a screening mass, it's existence would 
strongly suggest that a mechanism for the restoration of unitarity at very small x 
exists already at the classical level. 

We will address here the question of a screening mass in the classical theory 
quantitatively by solving stochastic difference equations on a two dimensional lattice. 
In Section 3, we will describe how we set up the problem and how one may use lattice 
perturbation theory to identify the weak coupling and strong coupling regimes of 
the theory. We define "reduced" correlation functions of gauge fields which are one 
dimensional projections of the original two dimensional fields. If a screening mass 
existed in the theory, these reduced correlation functions may be expected to have a 
very characteristic exponential fall off at large distances. 

For simplicity, we will consider an SU(2) gauge theory in our numerical work. We 
use the conjugate gradient method to solve the difference equations on the lattice. 
Details of the numerical procedure are also discussed in Section 3. In Section 4, we 
describe lattice results for the reduced correlation functions and compare them to 
the results expected from lattice perturbation theory in weak coupling. It is observed 
that in the weak coupling region, the numerical results reproduce to high accuracy 
the results of lattice perturbation theory. However, as one approaches the strong 
coupling region on the lattice, the number of the stochastic difference equations to 
which solutions can be found decreases and eventually, in the strong coupling region, 
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no solutions of the lattice equations exist at the desired O(o^) accuracy. In Section 
5, we will interpret these results and state our conclusions. 



2 Parton distributions as correlation functions of a 2— D 
field theory 

In the model of McLerran and Venugopalan, the partition function which describes 
the ground state properties of wee partons with a; << A^^^^ and transverse momenta 
qt« fm-i, is [| 

Z = J [dAtdA+][dtpU'tp][dp] 

exp(^iS + ig J d'^xA+{x)5{x~)p{x) - ^ J d^xtp'^{0,xt)j . (2) 

In the above, p is the valence quark color charge density. Also, the parameter p? ~ 
^1/3 f]^-2 jg ^YiQ average valence quark color charge squared per unit area. Since 
p^ is the only scale in the partition function above, the coupling constant will run 
as a function of this scale j^. If >> Aq^^j^, as will be true for very large nuclei, 
asip"^) << 1 and weak coupling methods can be used. 

If we integrate over the p fields first, we obtain an effective action for the wee 
partons with non-local propagators and vertices. Instead, the procedure followed in 
Ref. was to perform the p integrals last. In that approach, one needs to calculate 
the saddle point solution of the action for each p configuration to determine the 
classical background field. Any physical observable, such as a correlation function, 
is then obtained by evaluating it for the saddle point solution and then averaging it 
over all possible p configurations. The saddle point solution is nothing else but the 
solution to the classical Yang-Mills equations 

D^F^^ = 5J^ (3) 

in the presence of the external source J'^ = 5'^^ p{xt)5{x^). It was shown in Ref. Q 
that the solution of these classical equations of motion is 

^+ = 



3 



A- = 

= e{x~)a\xt) 



(4) 



The transverse components A^ , where j = 1, 2 further satisfy the equations, F12 = 
and V • a = gp, where the latter equation fohows from the Gauss' law. Since the 
fields A^ are thus gauge transforms of vacuum configurations with a gauge condition 
determined by the /^-configuration, one can write aj{xt) = {xt)V jU^xt) , where 
[/ is a unitary SU{3) matrix for QCD and an SU{N) matrix for a theory with 
colors. Substituting for aj in the gauge condition results in the stochastic differential 
equation 

V • UVU^ = -ig^p{xt) . (5) 

Using the solutions of the equation above, which are the saddle point configu- 
ration of the partition function in Eq. |2|, one can show that the classical correlation 
functions may be expressed (in matrix representation) as correlation functions of a 
two dimensional Euclidean field theory: 

{af{xt)af'm, = j^J[dp](u{xt)VUHxt)y^^(umuH0)y^^ 

X exp (^-^ J (fxtp''{xt)p''{xt)^ /I, (6) 

where the Gaussian random measure 

/ = j [dp] exp {-^ j ci2^i/(xi)/(x0) , (7) 

is all that is left from the original partition function. Note that the charges are highly 
localized in the transverse plane: 

{p%xt)) = ; {p%xt)p\yt))p = pH-'>5^^\xt - yt) . (8) 

In order to ensure that the valence quark color charge is confined to the transverse 
radius of the nucleus, we require that 

I d^xtp''ixt) = 0. (9) 
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This constraint was not stated in Ref. |^. In the momentum space this condition 
decrees that p"'{kt = 0) = 0. Thus the /cf = mode is excluded explicitly. 

To compute the correlation function in Eq. ^, we need to solve Eq. |5| to determine 
U = U{p) for each p configuration. We were unable to find an analytic solution 
to this highly non-linear equation for all values of the coupling. In the following 
and subsequent sections we will discuss the analytic weak coupling solution and the 
numerical solutions of this equation on a two dimensional lattice. 

For completeness, let us recall that the relation between distribution functions 
and the correlation functions above is straightforward and is discussed explicitly in 
Ret. §: 

where the trace is over both Lorentz and color indices. 

It was argued in Ref. |^ that the distribution function has the general form 

1 _(iV?-l)li^,,,/„.^.,_ 



7ri?^ dxfPkt vr^ x as 

where H{k'l / a^^p?) is a non-trivial function obtained by explicitly solving Eq. ^ 
The effective coupling constant of this theory was believed to be asfJ-/kt and that 
in the "weak coupling" limit asfJ- << h, H{k'l / Ogp?) a|/i^//cf, recovering the 
Weizsacker-Williams result scaled by p? . It was also conjectured that the function 
H would have the form Ogp? /{k'l + M^), where M = casn is a screening mass which 
is a constant c times the dimensionful scale asp- 

Interestingly, the problem formulated above is analogous to the problem of the 
critical behaviour of Ising~like models coupled to a random magnetic field. As dis- 
cussed by Parisi and Sourlas the partition function in that case has a structure 
identical to Eq. |2| albeit they only discussed the case of a scalar theory. It was argued 
in Ref. |15| that the singular behavior of the theory near the critical point was best 
described by correlation functions analogous to those in Eq. |^. The remarkable result 
of Parisi and Sourlas was that their scalar version of Eq. ^ could be written as correla- 
tion functions of a theory which is identical to the original theory without the random 
magnetic fields but va. D — 2 dimensions, where D is the dimensionality of the original 
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theory. This dimensional reduction is a consequence of a hidden supersymmetry of 
the expression analogous to Eq. The gauge theory analogue of this symmetry is 
nothing other than the well known BRST symmetry. Unfortunately, Parisi-Sourlas 
dimensional reduction will not apply to Eq. |^ because the analog of their scalar field 
is the compact field U as opposed to the gauge field a. 

Therefore, in order to test our conjecture about the existence of a screening mass 
in the the strong coupling domain, we address the question numerically by formulating 
the problem on a lattice. This also enables us to define the weak coupling limit more 
precisely. As we will discuss in the following section, the above statements about 
weak coupling are modified somewhat by the precise formulation of the problem on 
the lattice. The effective coupling of the theory is indeed asu/h but only for discrete 
multiples of kt = 2tt/L. Here L is the transverse size of the nucleus. Correspondingly, 
one obtains a discrete version of the Weizsacker-Williams result for weak coupling 
{g^ liL << 5 as we will show) by using lattice perturbation theory. However, since the 
limit L ^ oo is synonymous with strong coupling, the Weizsacker-Williams result of 
Ref. for continuous transverse momenta will not be recovered. 

3 The 2-D Theory on the lattice 

As discussed in the previous section, to compute correlation functions of the two 
dimensional field theory, we need to solve stochastic differential equations Eq. |5| for 
an arbitrarily large coupling. We intend to do this numerically by introducing a 
spatial lattice. The lattice spacing a serves as an ultra-violet regulator. Indeed, 
without such a regularization the functional integrals in Eq. |6| are not well defined 
since the correlations of the p-fields are proportional to a (5-function. Introducing 
the lattice, one sees from Eq. |^ that each p"(a;) is /i/a times a Gaussian random 
number of unit variance. Approximating the circular transverse side of the nucleus 
of diameter L by a square of length L, one sees it to be a x grid of lattice points 
with L = Na. The continuum limit consists of taking o ^ and N ^ oo such that 
L is held constant. The further removal of the infra-red regulator can be achieved 
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by taking L ^ oo, which was the hmit in which the weak couphng considerations of 
Refs. [H and Q led to their computational scheme of the low-x parton distributions. 

Using the unitarity condition on the [/-matrices, the stochastic differential equa- 
tion in Eq. ^ can be re-written as 

(UV^U^ -V^U -U^) = -2ig'^p, (12) 



On the lattice, finite differences replace the derivatives: 



„, M{xt + aj) + UHxt-aj)-2W{xt) 
V2[/t=^A _ L + o{a^). (13) 



i=l,2 



The labels j = 1,2 refer to the orthonormal directions on the lattice and aj denotes 
a displacement by a single site, i.e., by distance a, in the j'th direction. The resultant 
stochastic difference equation form (to O(a^) accuracy) of Eq. ^ is 



U{xt) J2 {uHxt + aj) + UHxt-aj: 

J=l,2 



h.c. + lig"^ fia p{xt) = . (14) 



In the equation above, h.c. deontes hermitean conjugate and we have scaled p 
fip/a. This has the advantage that the Gaussian random measure defined in Eq. |^ 
is now independent of the lattice spacing a and the dimensionful parameter p. It is 
redefined to be 

exp(^-i5:/(:r,)/(x,)^ . (15) 

The rescaled p on the lattice satisfy the following equations in analogy with their 
continuum version: 

ip'^ixt)) = ; {p'^{xt)p\yt))p = ^'^^d • (16) 

The zero net color charge constraint naturally becomes a sum on the lattice and 
p'^{kt = 0) = is true on the lattice as well. The only coupling this lattice theory 
has is the dimensionless g^pa and the scale for the theory is provided by the nuclear 
transverse size L. Physical quantities can therefore be obtained as a function of g^pa 
or equivalently g^pL. 
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In computing correlation functions on the lattice, we will find it most convenient 
to study correlations of one dimensional projections of the two dimensional gauge 
fields. These "reduced" gauge fields are defined as 

aj{X2) = • (17) 

This one dimensional projection sets the momentum ki = 0. If there exists a mass 
gap M in the theory, then 

(a;;(x2)a;;(y2)) = Aexp(-M |x2-y2|) , (18) 

for sufficiently large \x2 — 2/2! (to avoid influence of excited states, if any). An expo- 
nential fall off of the correlations of these reduced gauge fields would therefore be an 
unambiguous signature of a mass gap in the theory. 

3.1 Weak Coupling Limit 

Since the above mentioned function and therefore the correlation functions we 
wish to obtain, have earlier been obtained in weak coupling limit of the 2-D theory, 
it will be instructive to first calculate them in the weak coupling limit on the lattice. 



For small enough g^fia, Eq. 14 clearly admits a solution for U{xt) which is close to the 
identity matrix for all xt modulo a global gauge rotation. Writing the matrices U in 
terms of the the generators of the SU {N) gauge group as U{xt) = explig"^ ^a-(j){xt)) , 
with cl){xt) = J2k^^{^t) ■ T^, one sees that weak coupling implies that g^/xa • << 1. 
One can therefore expand the field U as 

U = l + ig^fia-cf>-^ig^fia)'' -^^ + ■■■ (19) 



Keeping terms of only lowest order, Eq. 14 becomes 



[<^{xt + Oj) + Hxt - aj) - 2(l){xt)] = p{xt) . (20) 

Thes equations can be solved by Fourier transforming the fields (j) and the sources p 
: Let 

1 ^ ^ / \ ~ /27rn\ 

^^^^) = 772 (2^^" • xt/L) 4> ( — I , (21) 



with n = (ni,n2) and — (A^ — l)/2 < nj < (A^ — l)/2 (assuming to be odd). 
Defining similarly the Fom'ier transform of pixt), we obtain the following solution of 
the Eq. M. 



\~) ~ ^=,,2 hs(^) - 1 



Substituting back in the equation for (j){xt), we have 



1 



2iV2 



E 



exp (27rm • Xt/L) 



ft Y.j=i: 



cos 



' 27rn, a % 



27rn\ 



In this leading order of weak coupling, aj{xt) = gfi [(/'(xt + aj) 



(22) 



[xt - a 



(23) 
and 



the one dimensional projections of the a fields, are easily computed to be 



al{x2) 
"2 (3^2) 



0. 

2N^ 



"2 



[cos 



27rn2a 

L j f2mn2X2\ 

exp 



1 



\ L 



(24) 



Here prime denotes the exclusion of the n2 = due to the total vanishing charge 
condition. To obtain the "reduced" correlators, we take the product of the fields 
and take the average over the p fields. Using the relation between p and p, and the 



Eq. 16, one can easily show that p satisfies similar equations as well except that its 
two point correlation function has an extra factor of N"^: 



Using the relation above and after some simple algebra, we obtain 



(25) 



{aUx)aUx')) 
{al{x)al{x')) 



0. 

y p \ " 

2Af2 ^ 



sin2 (2z[£2a 



In the continuum limit of a — > and N 



n2=l [cOs(^ 
CXD, 



cosl '""^^^-"^^ ). (26) 



(a2«2)a- 



2^2 



1 "-2 

712 = 1 ^ 

By constructing similar "reduced" correlators for the calculations of Ref. , one can 
easily see that our results are very similar to theirs, except that our expression above 



(27) 
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still has a finite L. Consequently, only discrete momenta are allowed in our sum and 
the lowest allowed momentum is 2tt/L. A naive L — > cxd yields identical results to 
those of Ref. [Q] but it turns out that this limit is not allowed. 

In order to see why it is so, it is necessary to go back to the weak coupling 
condition g^^a ■ (j) « 1- Using the solution for <j){xt) obtained above, one can 
translate this condition into 

1 ^ exp (27rin • x*j/L) 



Taking further the continuum limit, one obtains 

One can now see that the L ^ oo limit will violate the above condition even if 
one ignores the possibly logarithmically divergent factor in the curly bracket in that 
limit. The weak coupling condition thus constrains L to stay finite and small. The 
expression in the curly brackets can be evaluated numerically and typically the largest 
values are ~ 0.2 if one keeps L finite. The weak coupling condition for a finite size L 
is then (approximately) 

g'^fiL«5. (30) 

The correlation function {a'[a'[) becomes non-zero in the next-to- leading order of 
the expansion, when 

a/. = gy^^cp + ig^ {<P^f.<P - (V^0)0) . (31) 

For the sake of brevity we have used here the continuum notation to denote the finite 
differences. Using this expression, the reduced correlator can be computed straight- 
forwardly. The final expression is fairly tedious (involving the Gaussian average of 
four p fields). The key result of this computation is that in the continuum limit 
a ^ 0, 

{alal) ^g'p\g'fiLf, (32) 
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i.e., the correlation function grows as {g^iiL)'^ . In the next section, we wiU compare 
the results of our lattice computation with these lattice perturbation theory results 
in the weak coupling region g^^L << 5. 

The weak coupling calculations above were done by introducing an ultra-violet 
cut-off, the lattice spacing a. Since the final result shows sensitivity only to the 
infra-red regulator, namely the size L, one can ask whether these results can be 
derived without introducing the lattice in x-space at all. The answer turns out to be 
affirmative. One can easily show that the continuum problem can be formulated in the 
momentum space of a finite square box of length L. Due to the discrete momentum 
spectra, the corresponding /5-measure is then well defined. One then solves Eq. ^ by 
first expanding U{x) and then Fourier transforming the resultant equation. The final 
results, of course, remain unchanged when compared with the a — > limit above. 

3.2 Numerical Method 



In order to obtain a result for the correlation functions of the a''-fields which is free 
of the infra-red cut-off, one has to take the limit L ^ oo. Since it thus necessarily 
takes one out of the weak coupling region, we now turn to the procedure we used 



to solve the stochastic difference equations in Eq. 14 numerically. To simplify our 
computations and as a test, we choose to work with the gauge group SU{2). No 
qualitative differences are anticipated with regard to the existence of the mass gap 
as a result of this simplification. Writing an SU(2) matrix U as agl + iT^au, we can 



write the first term of Eq. 14 as 



i=i,2 



6oI + iT%k ■ 



(33) 



Here / is the unit 2x2 matrix and T^,k = 1, 2, 3 are the Pauli matrices. The coeffi- 
cients Ofc satisfy the unitarity condition, J2k=o '^k~^ coefficients do not. 



Eq. 14 can now be re-expressed as bk{xt, xt ± aj) +g'^fia ■ Pk{xt) = 0. In order to solve 
these coupled nonlinear equations, we minimize the function F, defined by 



F = Y^I^Y(bk{xt,Xt±aj) + g^fia ■ pk{xtfj +(^al{xt)-lj |. 



(34) 
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Minimizing F is equivalent to solving Eq. 14 for each lattice point and color charge 
(3 N"^ equations) while simultaneously imposing the unitarity condition U W = 1 at 
each point on the lattice. The latter is done by the second set of terms in F. Note 
that F is a sum of squares of real numbers with zero as its possible absolute and 
desired minimum. A lack of solution will be signalled by large values of -Fmin for the 
absolute minimum. 

We use a multi-dimensional conjugate gradient method, described in the subrou- 



tine FRPRMN and its associated subprograms in Ref. |14], to minimize F to an 
accuracy better than O(o^) as dictated by accuracy of the original lattice equations. 
The zero net charge condition compels us to use periodic boundary conditions in 
accordance with the Gauss' law. We investigated both ordered and random starts 
for the initial guesses for the f/'s. Each iteration consisted of choosing the source 
distributions randomly over the entire lattice in the momentum space such that 1) 
/5(0, 0) = 0, 2) p*{k) = p{—k) and 3) both the real and imaginary parts of each p{k) 
were random Gaussian numbers with variance l/\/2- The p-distribution was then ob- 
tained by an explicit Fourier transformation. Using the conjugate gradient method, 
the set of C/'s for the absolute minimum was found. If the minimum was O(a^) or 
smaller, then the matrices U = U{p) were used to compute the correlation functions 
on the lattice. We obtain from the relation 

Tk-a^^ = -^Im (u{xt)U\xt + a^)) . (35) 

We have also checked that the symmetric difference definition for the derivative yields 



the same result. Just as in Eq. 17, we define the "reduced" gauge fields and 
compute correlators by taking the product of these gauge fields. This procedure 
is repeated over several iterations, typically a few hundred, and the the correlation 
function is averaged over these iterations. The errors are determined in the usual 
way by computing the standard deviations. Note that the sets of p's in successive 
iterations are totally independent and one thus has negligibly small auto-correlations. 
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4 Results 



In view of the facts that that the couphng for the lattice theory above is g^^ia, and 
that none of the three quantities in this expression occur independently, we chose 
to set g^/i = 1 in our simulations and varied g^^a by varying the lattice spacing 
a. Simulations were performed for a range of lattice sizes, ranging from = 21 to 
N = 211, and for values of g^^L ranging from 0.5 to 20. Typically 200 iterations 
were performed, each consisting of an independent set of the /^-distributions, unless 
stated otherwise. Noting from Eq. ^ that sets the scale of the correlation 

functions, and using the definition in Eq. one can show that the factor {g'^fj.a)'^ 
relates the dimensionless lattice correlation function to the physical a^-correlations. 
We therefore show the results for the latter in the units of ^^/i^. 

In Fig. la we show the results of our computation for {a^j{x)a^j{0)) /g'^n'^, j = 1, 2, 
plotted as a function of the dimensionless distance x/L for a small value of g^fiL = 0.5. 
The lattice size was 41 x 41. Also shown are the analytic weak coupling results of 



Eq. 26 for this lattice size which agree rather well with the direct computation . One 
also sees clearly that the aj^-correlation is very small compared to the Q2-correlation. 
In fact, the former is consistent with zero on the scale of this plot. This, together 
with the excellent agreement with the weak coupling result, reassures us that 1) for 
small g^i-iL the assumptions made in deriving the weak coupling results are indeed 
justified and 2) our numerical procedure works fine. 

Fig. lb further shows that these results are indeed the continuum results. It 
displays the results for the same g^fiL but on = 21 and 41 lattices. The results for 
the j = 2 correlation are again displayed as a function of the dimensionless distance 
x/L and the results are seen to be lattice size independent. One may wonder why the 
correlation function is negative, given that it can now be thought of as a continuum 
property. The weak coupling result of Eq. ^ provides a hint for understanding this. 
The leading term in it is negative for x — x' = L/2 and the successive terms alternate 
in sign and become progressively smaller in magnitude. Thus for any finite L, the 
correlation function will be negative midway if one is in the weak coupling domain. 
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Although the {ai{x)ai{0)) correlation function appears to be zero at all x/L in 
Figs, la and lb, it has an interesting structure as well. As Fig. 2 shows for = 21 
and 41 lattices, this correlation function decreases monotonically as x increases but 
remains positive all through. We will later compare this behavior for larger g^fiL but 
it is interesting to note this difference with the leading order weak coupling behavior. 
As remarked earlier in Section 3.1, we do expect a non- vanishing contribution to 
it from the next-to-leading order contribution and a detailed examination of it also 
reveals it to be positive definite. 

Having tested both the weak coupling limit and the conjugate gradient method 
on the lattice, we increased the g^^L, first by retaining the same lattice size of A = 21 
and then increasing it as well up to A = 71 such that the lattice spacing stayed at 
a ~ 0.1. This value was determined by making runs on the A = 21 lattice for various 
a and by checking that the errors due to finite a remained small. For the rest of 
our numerical work we have attempted to stay close to this value of a; increasing 
thus the lattice size A^ in order to increase g^fiL. Fig. 3 displays the results for the 
a2-correlation functions in the units of g'^fj,'^ as a function oi x/L for g'^fiL = 0.5,1, 
2, 3, 4, 5, 6, and 7. It appears to remain almost independent of g^^L until it reaches 
our estimated region of the validity of the weak couping theory: for g^^L > 5 the 
correlation function tends to be less and less negative as g'^fiL increases. This signals 
a departure from the leading order weak coupling result which, as seen in Eq. is 
independent of g^fiL when viewed as a function of the dimensionless variable x/L. 



An obvious source of the departure from Eq. |2g are higher order contributions. If 
these tedious terms are indeed responsible for it then one expects a growth in the 
aj-correlation as we argued in Section 3.1. 

Fig. 4 exhibits the a^^-correlation function in the units of g^fi^L'^ as a function of 
x/L for g^fj.L up to 5. They do indeed group together to suggest a universal curve, 
and thus confirm the rise of this correlation function as {g'^fiL)'^. Fig. 5 demonstrates 
this in another way and also suggests that g'^fJ.L ~ 5 is the boundary of the weak 
coupling region. What is shown there is the x = value for this correlation function 
in the units of g'^iJ? as a function of g'^iiL, both before and after scaling out the factor 
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{g'^fiL)'^. Note the scale of both the axes. A hnear rising curve is thus an indication 
of the power law which seems to be consistent with the power two. What can also 
be inferred from this figure is a small trend to push this power up as one goes above 
g^fiL ~ 5. A priori, such a behavior could also be due to yet more higher order terms. 
However, these results also suffer from a further defect. 

For the larger values of g^fiL, one sees the Fynm slowly creep up and go beyond 
the O(a^) level. Indeed, typically one fails to obtain any acceptable minimum at that 
level for about 10-15 % of the iterations. This should be contrasted with the small 
g^fJ,L case where -Fmin was a lot smaller than O(a^) for each iteration. Increasing the 
g^fiL even further, this becomes worse very quickly and by g^^L = 10 no minimum 
exists at that accuracy. 

In order to better understand the reason behind this, we show in Fig. 6 a normal- 
ized histogram plot for g^fiL = 6, 7, 10 and 20. These runs were made on = 61, 
71, 111 and 211 lattices and the latter two have very few iterations, being 11 and 6 
respectively. All the corresponding minima were too high compared to O(a^). The 
parameter R in Fig. 6 is defined as follows. Defining R^{xt), k=l, 2 and 3, to be the 
terms on the LHS of Eq. 14 divided by a^, one sees that R^{xt) = 0, for all xt and 
k, is the desired solution. A value of R^{xt) 7^ measures how far away from the 
desired minimum (found by minimizing F) is the solution for that value of xt and k. 
What Fig. 6 depicts, for different values of g^fJ,L, is the fraction of the 3A^^ equations 
which have, for the best minimum of F, R given by the value on the x-axis. We have 
checked that the similar histogram plots for the weak coupling region have only the 
bin near zero occupied, i.e. they peak sharply at zero. What one sees in Fig. 6 though, 
are increasing deviations away from R = ~- fewer and fewer of the 3A^^ equations 
are being satisfied at the required level of accuracy. Noting that R = 1/a, which is 
~ 10 for these runs, corresponds to the equations not being satisfied at 0(a'^) level, 
one finds that the minimum of F has increasingly many equations like that. This is 
thus an indication that for g^fj.L > 10, no solutions to the stochastic equations exist 
at 0(0^). 

Our results therefore suggest the following: when g^fxL < 5, the correlation 
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functions computed directly agree very well with the expectations from lattice per- 
turbation theory. For intermediate values, g^^L ~ 5, the lattice results still agree 
reasonably well with the analytical lattice expressions but one notices an increasing 
trend of lack of solutions to more and more equations at the level. For larger 
values of g^^L > 10, no solutions exist at that level. 

The absence of solutions as we increase g^^L was unexpected. Increasing g^^L 
for a fixed value of a is equivalent to merely increasing the number of sites A^. In 
other words, the number of equations has been increased but the structure of the 
equations and the coupling is unchanged. Why then are there no solutions as we go 
beyond g'^^L ~ 5? One way to understand this is as follows: since U = ex.p{ig'^fia ■ 
(j)) and in weak coupling cp oc N, increasing g^fiL will cause the U matrices to 
deviate increasingly away from identity. However, since g^fJ-a is unchanged and it 



remains small, and since p remains 0(1), Eq. 14 will still prefer the U to be close 
to identity. The ensuing mismatch will thus result in a lack of solutions for large 
g^fiL. It is possible that our lack of solutions is because our boundary conditions are 
too restrictive. However, because we need to satisfy Gauss' law in two dimensions, 
periodic boundary conditions appear to be the appropriate physical choice. 

The absence of solutions for g'^fiL > 5 is a serious problem for the classical 
theory discussed in Ref. ||5|. Not only because the conjectured scenario of a screening 
mass needed the coupling to be strong but also because the removal of the infra-red 
regulator pushes one in that region. Considering that even the very large nuclei will be 
finite in size, one could check whether the condition above is physically acceptable. If 
we take L ~ 2 A^^'^, then the weak coupling condition g^fiL < 5, holds only for very 
small A. This is due to the fact that /x ~ A^^^ and the coupling g^ is also evaluated 
at the scale /i. Thus this condition contradicts the weak coupling assumption of the 
four dimensional theory in Ref. ^ which is expected to be valid only for very large 
A. 
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5 Summary and Outlook 



In Ref . 0] , a QCD based model was formulated to study the properties of low x, wee 
partons in large nuclei. For very large nuclei, it was argued that the problem could 
be formulated as a weak coupling, many body problem. In Ref. |^, it was shown 
that the classical saddle point solution of the model could be expressed as a two 
dimensional Euclidean field theory with the dimensionful coupling 05/^. Computing 
correlation functions in this 2-D theory required the solution of highly non-linear 
stochastic differential equations in the presence of a Gaussian random source. No 
analytic solution of these equations was found in Ref. 1^. 

However, it was argued that classical distribution functions had a Weizsacker- 
Williams distribution at large momenta kt » asf^- It was also conjectured that 
at smaller momenta, in "strong coupling", the theory acquired a screening mass 
M ~ ctsf^, which regulated the growth of the distribution function at small kt- The 
existence of a screening mass would be suggestive of a weak coupling, albeit non- 
perturbative, restoration of unitarity already at the classical level. 

In this paper, we have discussed the analytic weak coupling and the numerical 
solutions of the stochastic differential equations on a two dimensional lattice. For our 
numerical work, we made the simplifying assumption of two colors and investigated 
an SU{2) gauge theory. With lattice perturbation theory as our guide, we identify 
g^^iL << 5 as the weak coupling condition. Our numerical results on the lattice agree 
very well with lattice perturbation theory for these values. For larger values of g'^^L, 
no solutions are found which satisfy the stochastic equations at the required level of 
accuracy. Thus not only is a screening mass absent but a further implication of this 
result is that the classical theory is ill defined in the infrared. Furthermore, if we 
identify L ~ 2 A^^^ fm, the lattice weak coupling condition is satisfied only for very 
small A. This limit appears to contradict the weak coupling limit in the full theory, 
which is expected to be valid only for very large nuclei. In sum, our work suggests 
that the classical theory in Ref. |||] is seriously flawed. 

Recently, McLerran and collaborators have proposed that the original clas- 
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sical theory is flawed because the authors in Ref. failed to properly solve the 
Yang-Mills equations for the transverse components A'' of the classical field. These 
are determined through the equation 

\/id+A' + A,x d+A' = gj+ . (36) 

Ref. ||5| argued for a solution of the form in Eq. ^: A^ = a'^9{x~). If one then 
ignores the commutator terms, because it involves fields at the same x^, one then 
obtains V • a = gp^ which gives us the stochastic equation we solved for U using 
ai = UViW /{—ig). The authors of Ref. argue that the cross product term 
above cannot be dropped because of its peculiar singular structure. They argue that 
the source term must be regularized so that instead of being a (5-function in x~ , 
the charge density p depends on the spacetime rapidity y = — log(x~). The above 
equation is then re-written as 

= gp{y,xt) , (37) 

ay 



where Di is the covariant derivative. It is claimed in Ref. [16| that this equation 
can be solved exactly and the correlation functions computed analytically. The 
distribution functions have the Weizsacker-Williams form for large transverse mo- 
menta, dN/d^kt ~ • For small transverse momenta, it has the logarithmic form 

dN/d'^h ~ \og{k^Jx{y,k^t))- Here xivM) = IZx^^^, dy' f^Hv' ,Q^)- We refer the 
reader to Ref. llql for the details of their calculation. 
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Figure Captions 



Figure 1: (a) The correlation functions T^{x) = {a^j{x)a'j{0)) as a function of x/L for 
j = 1 (circles) and j = 2 (crosses) for g^fiL = L = 0.5 and L = Ala. The continuous 
line is the weak coupling result of Eq. (b) Same as Fig. la but for L = 21a and 
L = Ala and for only. 

Figure 2: The Oj^-correlation function as a function of x/L for g^fiL = 0.5 and = 21 
and 41 lattices. 

Figure 3: The a2-correlation function as a function of x/L for g^^L = 0.5 and 1, 2, 
3, 4, 5, 6, and 7. The lattice sizes can be found in the text. 

Figure 4: The a^^-correlation function in the units of g^fi'^L^ as a function ol x/L for 
g^fiL = 0.5 and 1, 2, 3, 4, and 5. The lattice sizes can be found in the text. 

Figure 5: The a^-correlation function at x = as a function of y^/xL without (crosses) 
and with (squares) a division by (g'^fiL)^). 

Figure 6: The histograms of the fraction of total equations solved at an accuracy 
* as a function of R. These are plotted for the following values of g^i^L: 6, 7, 
10 and 20. For details, see the text. 
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